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A simple geometrical model is presented for the gravity-driven motion of a single particle on a 
rough inclined surface. Adopting a simple restitution law for the collisions between the particle and 
the surface, we arrive at a model in which the dynamics is described by a one-dimensional map. 
This map is studied in detail and it is shown to exhibit several dynamical regimes (steady state, 
chaotic behavior, and accelerated motion) as the model parameters vary. A phase diagram showing 
the corresponding domain of existence for these regimes is presented. The model is also found to be 
in good qualitative agreement with recent experiments on a ball moving on a rough inclined line. 
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I. INTRODUCTION 

Several experimental studies have recently been 

conducted on the problem of a single ball falling under 
gravity on a surface of controlled roughness. These works 
have revealed interesting new aspects of granular dynam- 
ics that are not yet fully understood. Three distinct dy- 
namical regimes have been identified fll-Ql as the tilting 
angle increases. For small inclinations there is (i) a decel- 
erated regime where the ball always stops, then comes (ii) 
an intermediate regime where the ball reaches a steady 
state with constant mean velocity, and for larger inclina- 
tions the ball enters (iii) a jumping regime. Computer 
simulations 1^-^ have confirmed these results, particu- 
larly those concerning regimes (i) and (ii). A theoretical 
model has also been proposed in which steady-state 
solutions (but no detailed dynamics) can be obtained an- 
alytically. More recently, a one-dimensional map H] has 
been introduced to study the jumping regime. This map 
in its simplest version is linear, and to obtain non-linear 
behavior one has to vary spatially the properties of the 
rough surface in which case the model become inac- 
cessible analytically. 

In this Paper we present a model for a single particle 
moving under the action of gravity on a rough surface 
of specified shape. Within this setting we will give a 
detailed analytical description of all possible dynamical 
regimes. Although the model we study is simplified, its 
predictions are in good qualitative agreement with the 
experimental findings. 

Roughly speaking, our conclusions are as follows. 
There is (i) a sharp transition (as the surface inclination 
increases) from a regime of bounded velocity to one of 
accelerated motion. Within the region of bounded veloc- 
ity various dynamical regimes are possible. First there is 
(ii) a range of inclinations for which the dynamics always 
has a unique attractor. For higher inclinations two other 
phases exist: (iii) a region where we have co-existing at- 
tractors for the dynamics and (iv) a region where instabil- 
ities give rise to chaotic behavior. For a fixed (sufficiently 



large) inclination a transition to the chaotic region will 
take place as the nature of the collisions between the par- 
ticle and the surface becomes highly inelastic. Although 
our results are derived here in the context of a simple 
collision rule, it can be shown |^ that they remain valid 
for a wide class of restitution laws. 

The paper is organized as follows. In Sec. II we de- 
scribe the model and study in detail its dynamical prop- 
erties. The main results of this Section are then sum- 
marized in the phase diagram shown in Fig. ^ In Sec. 
Ill we carry out a comparison between the model pre- 
dictions and the experimental findings. In particular, we 
argue that the jumping regime seen in the experiments 
might correspond to a true chaotic motion, as predicted 
by the model. Finally, in Sec. IV we collect our main 
conclusions and present further discussions. 

II. THE MODEL 

In our model, which is shown in Fig. 0, the rough sur- 
face is considered to have a simple staircase shape whose 
steps have height a and length b. For convenience, we 
choose a system of coordinates in such a way that the 
step plateaus are aligned with the x axis and the direc- 
tion of the acceleration of gravity g makes an angle (p 
with the y axis. A grain is then imagined to be launched 
on the top of the 'staircase' with a given initial velocity. 
In what follows, we will be concerned with the problem of 
a point particle falling down this 'staircase' and will thus 
not take into account any effect due to the finite size 
of the grain. Upon reaching the end of a step plateau, 
the particle will undergo a ballistic flight until it collides 
with another plateau located a certain number n of steps 
below the departure step (e.g., n = 3 in Fig. 1). Accord- 
ingly, we will refer to the integer n as the jump number 
associated with this flight. 

We will assume, for simplicity, that the momentum 
loss due to collisions is determined by two coefficients 
of restitution ct and e„ , corresponding to the tangential 
and normal directions, respectively. More precisely, if 
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V = {vx,Vy) denote, respectively, the components of the 
particle velocity parallel and perpendicular to the surface 
before a collision, then we will take the velocity v' = 
{v'^,v'y) after the coUision to be given by 



(la) 
(lb) 



where < et < 1 and < e„ < 1. 




FIG. 1. Model for a single particle moving under gravity 
on an rough inclined surface. 

In the present paper we will for simplicity discuss only 
the case e„ — 0; the advantage being that the model 
can then be described by a one-dimensional map. When 
e„ > the dynamics is governed by a three-dimensional 
map, the analysis of which is more complicated and will 
be left for forthcoming pubhcations [ p^ . 

We now derive the equations governing the dynamics 
of the model presented above. Let us denote by E the 
kinetic energy of the particle at the moment of departure 
for a given flight. We write E = ^mV^, where m is the 
particle mass and V is the launching velocity at the start 
of the flight (see Fig. |^). After this flight the particle 
will first collide with a step below, then slide along this 
step (recall e„ = 0), and finally take off again on another 
flight with initial kinetic energy E'. We suppose that the 
main energy loss is due to collisions and so we neglect 
the energy dissipation as the particle slides along a step, 
where it then moves with a constant acceleration g sin (p. 
Using simple arguments of energy conservation together 
with the collision conditions (1) and (2), one can write 
E' in terms of E. The result is 



E' = —me^v'^ + mg sin {nb — x), 



(2) 



where n is the corresponding jump number for the flight 
and X is the cc-coordinate of the landing point. It takes 
a simple algebra to show that at the landing point (x, y) 
we have the following identities: 
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(3a) 



.9 cos 2 

— - — 1 = na, 



g sin (f) T - 
g cos (j) T, 



(3b) 

(3c) 
(3d) 



where T is the flight time. 

It is convenient to introduce a dimensionless energy- 
like variable: 



£ = 



E 



mga cos <j> 



(4) 



Eliminating T from (|^) and inserting the result into (||), 
we obtain that the dynamics of the model in terms of the 
variable £ is given by the following map: 



£' = f{£,n) = n 



(^y/£jVi + t^^ + t i^T - t - 2y/Ffnj 

(5) 

where we have for conciseness introduced the notation 



tan (I 
: h/a. 



(6) 
(7) 



The parameter r above can be viewed as a measure of 
the surface roughness, with — corresponding to a 
perfectly smooth surface. As for the inclination parame- 
ter i, we need to consider only the interval < i < t for 
which non-trivial motion occurs. (Clearly, for t < the 
particle will always come to a rest, whereas for t > t the 
particle undergoes a free fall without ever colliding again 
with the ramp.) 

The flight jump number n appearing in Eq. (^) is de- 
termined from the energy £ according to the following 
condition: n is equal to the smallest integer such that 
n& — 2; > or, alternatively. 



n(T - t) 



i£ > 0. 



This means that £ falls within the interval /„ 



£ e ln{t) 



i(n-l)(r-t)2,ln(T-i)2 



(8) 



(9) 



Thus the function f{£,n), as defined by Eqs. (||) and 
(^), exhibits jump discontinuities at energy values £ = 
jfi (r — t)^, but each of its branches is smooth. This is 
illustrated in Fig. |^, where we graph the function (^) for 
et = 0.7, r = 3.7, and several values of the inclination t, 
For later use, we note here that the average velocity V 
between two consecutive fiights is given by 



V 



7lL 



T + 2E' /m~etv^)/g sin (j)' 



(10) 



where L = y^a^ + P and the second term in the denom- 
inator corresponds to the time during which the particle 
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moves on the ramp (see Fig. 1). If we now introduce a 
dimensionless mean velocity 



V 



\/ ag cos (j) 
then Eq. (|o|) becomes 

V^(1 + t2)/2 



V 



(1 - et)t + \/£V^ - et^t/^ 



(11) 



(12) 



In order to study the dynamical properties of the map 
above, we must first investigate the existence of fixed 
points. If wc denote by a fixed point with a jump 
number n, then will be the solution to the equation 



£n = f{£n,n). 



(13) 



In view of the homogeneity of the function f{£,n) [see 
Eq. (^] we write 



£n = n[zo{t)Y 



(14) 



where the quantity zo(t) no longer bears any dependence 
on n. Using Eqs. (^) and (|lj), Eq. ( p^ ) becomes 



(zo + t)2 = e?(<- 
whose positive solution is 



zo{t) 



-t + 



Tt 



(15) 



(16) 



Now a fixed point £„, as given in Eqs. ( |14| ) and (16), 
will exist if and only if £„ e In{t), where the interval 
Init) is defined in (|). Thus, as t increases, a fixed point 
with jump number n will be born when £„ equals the left 
endpoint of /„. Comparing Eqs. (^, (|lj) and (|l6|), we 
see that this happens at an inclination f„ such that 



zoiU. 



-tn + 



Tt, 



(r-t„). (17) 



This equation is quadratic in ^/t^ and can thus be easily 
solved. However, we shall not bother to give the result 
here and will simply mention a few important facts that 
follow from Eq. ([l7|). First, we note that ii = so that 
a fixed point with jump number n = 1 is always born at 
t = 0. Then, as t increases, fixed points with successively 
higher n will appear in an increasing sequence of inclina- 
tions {irij^i- Finally, we have that for t > too, where 
too — hm„_^oo ^ri, all fixed points cease to exist. Setting 
n = oo in Eq. ( [l7| ) we obtain for the limit point too- 



tryn — T 



l + et' 



(18) 



The appearance of this sequence of fixed points can per- 
haps be best visualized by referring to Fig. 0, where we 



plot the function f{£,n) at increasing values of t, with 
et and r kept fixed. For small t (lower-most curve in Fig. 
^ there is only one intersection with the 45° line, corre- 
sponding to the fixed point with n — 1. As t increases 
fixed points with successively higher n appear (second 
curve from the bottom). At t = too there are infinitely 
many such fixed points (second curve from the top) and 
after this all of them cease to exist (uppermost curve). 



f(£) 




FIG. 2. One-dimensional map f{£,n) for et — 0.7, r — 3.7 
and t = 0.2, 0.5, 0.7, 0.9 (from the bottom up). 

One can also show that for t > too we always have 
f{£,n) > £, whereas for < t < too there exists an en- 
ergy £* such that f{£,n) < £ for £ > £* (see, e.g.. Fig. 
§) . We thus conclude that for t > too the particle velocity 
will become unbounded for any initial condition, whereas 
for < t < too the velocity remains always bounded. 
In other words, at the critical inclination t = too there 
is a sharp transition (independent of initial conditions) 
from a regime of bounded velocity to accelerated mo- 
tion. In the region of bounded velocity, several dynam- 
ical regimes are possible, depending on the stability of 
the fixed points, as discussed below. 

The stability of a fixed point is determined by the 
parameter A = f'{£m n), where the prime denotes deriva- 
tive with respect to £, so that if |A| < 1 (|A| > 1) the 
fixed point is stable (unstable) ||ll|. Using Eqs. (|), (|l|) 
and (|l^), we obtain for the derivative A at the fixed point: 



A(t) 



1 



1 



l-^il^ei)t/T 



(19) 



Notice that A does not depend on n, thus implying that 
all existing fixed points £„ (for given values of the model 
parameters) have the same stability properties. More- 
over, since A is always smaller than unity, we see that 
instability can occur only if A(t) < —1. Let us then 
denote by tinst the inclination such that A(tinst) = — 1- 
FromEq. (|l|) we obtain that 



4(1 -e?) 



(20) 
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Thus the fixed points are stable for t < t-mst and unstable 

for t > tinsf 

If the fixed points are stable, the dynamics of the map 
will in general be attracted to one of the existing fixed 
points. For example, in the region of parameters such 
that < t < i2 < ^inst the particle will almost always 
reach a periodic motion where the particle falls by one 
step at a time, since in this case only the fixed point with 
n = 1 exists and is stable On the other hand, for 

t2 < t < tinst there are co-existing stable fixed points, in 
which case the final state (i.e., the fixed points to which 
the dynamics is attracted) will depend on the initial con- 
dition. Once the system has reached a given fixed point 
£n the particle will accordingly be moving with a con- 
stant mean velocity Vn whose value can be readily ob- 
tained by inserting Eqs. ( |l^ ) and into Eq. ([l^): 



When the fixed points are unstable (iinst < t < too), 
the particle motion becomes very irregular and no sta- 
tionary (periodic) regime is ever reached. This is illus- 
trated in Fig. ^, where we plot the jump number n as 
a function of time (iteration step) for two orbits in the 
region where the fixed points are unstable. In this figure 
we clearly see that the jump number fluctuates erratically 
around a mean value. We have computed the Lyapunov 
exponent for several values of parameters in the region 
of unstable fixed points and have found it to be positive 
for all cases studied, thus indicating that the motion is 
indeed chaotic in this region. 
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FIG. 3. The jump number n as a function of time (mea- 
sured in iteration steps) in the chaotic regime. Here et = 0.35, 
T — 3.73, and t = 1.65 (lower orbit), 1.74 (upper orbit). 

The different dynamical regimes displayed by the 
model above can be conveniently summarized in terms 
of a "phase diagram" in the parameter space (e^, t/r), as 
shown in Fig. ^. In this figure we plot the curves corre- 
sponding to too (solid line) and ti^st (dashed line) given 
by Eqs. ( p^ and (^, respectively. Also plotted is the 



curve representing the inclination t2 (dot-dashed line) at 
which the fixed point with n = 2 first appears. Thus 
in terms of the existence/stability of the fixed points 
the model displays the following four regions: (i) for 
< t < min(t2, ^inst) there is a unique stable fixed point, 
namely, that with n = 1; (ii) for t2 < t < min(tinst, ^oo) 
there are multiple stable fixed points (at least those with 
ri = 1 and n = 2); (iii) for t-mst < t < too all existing fixed 
point are unstable and chaotic motion is observed; (iv) 
for t > too no fixed point exists and the motion becomes 
accelerated. 

Another interesting feature in Fig. ^ is the fact that 
the chaotic regime appears when the collisions are highly 
inelastic (i.e., small et). In particular, for et > V2— 1 (at 
which point tinst equals too) the fixed points remain stable 
over their entire domain of existence. (The results shown 
in Fig. ^ are qualitatively different from the behavior seen 
in the model studied in Ref. where chaotic motion 
appears as the restitution coefficient increases.) 




FIG. 4. Phase diagram for the model. The solid line cor- 
responds to too, the dashed line to tinst, and the dot-dashed 
line to t2. 

III. COMPARISON WITH EXPERIMENTS 

In this section we wish to compare our model with re- 
cent experimental studies of a single ball moving under 
gravity on a rough inclined surface. In these experiments, 
first performed by Jan et al. jl] and later expanded by 
Ristow et al. , a rough surface was constructed by glu- 
ing steel spheres of radius r on a L-shaped flume. An- 
other steel sphere of radius R was then launched with 
a small initial velocity and its subsequent motion ana- 
lyzed. As the surface inclinations increases, the follow- 
ing three regimes are observed [Q : for small inclinations 
the bead always stops (regime A), then comes a range 
of inclinations for which the ball reaches a steady state 
with constant mean velocity (regime B), and beyond this 
point the ball starts to jump (regime C). In Fig. || we 
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show data taken from Ref. |Q] for the ball mean velocity 
as a function of sm9, where 6 is the inclination angle 
with respect to the horizontal direction. As discussed in 
Ref. Q, the change in trend observed in the data as 
increases (for a given value of R/r) marks the beginning 
of the jumping regime. 

The regime B seen in the experiments corresponds in 
our model to a stable fixed point with n — 1, for in this 
case the particle reaches a periodic motion where it falls 
one step at a time (as in the experiments). In order to 
compare our model more closely with the experiments let 
us first express the mean velocity Vi (at the fixed point 
n — 1) in terms of the angle 9, where 6 = (f> + Tr/2 — a 
(see Fig. 1). Setting n = 1 in Eq. (|2l|), returning to 
dimensionful units via Eq. (|lTl), and expressing the final 
result in terms of 9, we obtain 



Lg{l + et) 



-.1/2 



I 2(1 -et) 



VI- 



sm y — r 



(22) 



(We remark parenthetically that a similar expression can 
be obtained heuristically if one introduces an effective 
sliding friction in addition to inelastic collisions; see Refs. 
[0,1). Our formula follows however from a pure collision 
model.) 

We have fitted the expression (^2|) to the experimental 
data shown in Fig. ^ — the corresponding results being 
displayed as solid curves in this figure. In our fitting 
procedure, we took L = 2r = 1 cm [Q, g = 980 cm/s^, 
and best-fitted the parameters r and et for each data set 
considering only points in regime B. As we see in Fig. |^, 
the model prediction for the dependence of V with 9 is in 
a good agreement with the experimental data (in regime 
B). _ 

The jumping regime observed in the experiments, on 
the other hand, would correspond in our model to the 
region of unstable fixed points, since in this case the par- 
ticle jumps erratically never reaching a steady state (see 
Fig. H). This analogy might then provide a possible ex- 
planation for the change in trend observed in the exper- 
imental data for large inclinations. To see this, consider 
the region of small et in the phase diagram shown in Fig. 
^. As the inclination t increases (for a given et) the sys- 
tem goes from a region of stable periodic motion (with 
n = 1) to a regime of chaotic jumps, in close resemblance 
to the experimental transition from steady-state to the 
jumping regime. 

To probe this analogy further, we illustrate in Fig. ^ 
the behavior predicted by the model for the mean veloc- 
ity y as a function of sin 9 in the region of small et . In 
this figure, the solid curve corresponds to the expression 
( p^ ) for Vi, up to the point where the fixed point goes 
unstable, and the crosses are computed values of V in 
the ensuing chaotic regime. Comparing Fig. || with Fig. 
^, we see that the change in behavior predicted by the 
model at the onset of instability is in qualitative agree- 



ment with what is observed in the experiments (for small 
values of R/r) as the ball enters the jumping regime. Of 
course, more detailed experiments are necessary to ver- 
ify whether chaotic motion does indeed take place in the 
jumping regime. 




FIG. 5. Mean velocity V (cm/s) as a function of sin^. 
Points are experimental data taken from Ref. [4] for R/r = 2 
{+), 1.5 (*), 1 (A), 0.8 (O). Solid curves are theoretical fits 
[Eq. (|2^)], ending near the last data point considered in the 
fit. Fitted parameters are (et,T) = (0.72, 33.18), (0.64, 21.09), 
(0.41, 10.17), (0.27, 7.07), from left to right. 
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FIG. 6. Same as figure ^ for our model with et =0.1 and 
T = 1. The solid curve corresponds to Eq. (^^ whereas the 
stars give the computed mean velocity in the chaotic regime. 



IV. CONCLUSIONS 

We have studied a simple geometrical model for the 
gravity-driven motion of a single particle on a rough 
inclined line. In our model the rough line was chosen 
to have a regular staircase shape and a simple collision 
law was adopted. With these simplifications the dynam- 
ics is described by a one-dimensional map that is quite 
amenable to analytical treatment. Summarizing our find- 
ings, we have seen that our model displays the following 
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four dynamical regimes: 

1. for < < < min(i2, iinst) there is a unique stable 
fixed point. 

2. for t2 < t < min(iinst, ^oo) the system has multiple 
stable fixed points. 

3. for tinst < t < too the fixed points are unstable and 
the dynamics is chaotic. 

4. for t > too no fixed point exists and the motion 
becomes accelerated. 

Here the parameter t measures the surface inclination 
and the quantities t2 , iinst , and too separating the differ- 
ent regimes are given in terms of the other two model- 
parameters, namely, the restitution coefficient et and the 
roughness parameter r. These regimes are indicated in 
the phase diagram shown in Fig. y. Furthermore, it can 
be shown that the above conclusions, which were de- 
rived in the context of a simple collision rule, remain valid 
for a wide class of tangential restitution laws. 

Despite its simplicity, our model does provide a theo- 
retical framework within which the generic behavior seen 
in experiments on a ball moving on a rough surface can 
be qualitatively understood. For example, the model 
successfully predicts the existence of several dynamical 
regimes that are also observed in the experiments. In 
particular, the predicted functional dependence of the 
mean velocity with the inclination angle 9 (in the steady- 
state regime) is in good agreement with the experiments. 
Moreover, the model provides a possible explanation for 
the change in trend seen in the experimental data as the 
ball enters the jumping regime. We have suggested that 
this jumping regime might correspond to a chaotic mo- 
tion, as so happens in the model. Clearly, more experi- 
mental studies are required to investigate this interesting 
possibility. 
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